Calculation of the hyperfine structure of the superheavy elements Z=119 and Z=120 
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The hyperfine structure constants of the lowest s and P1/2 states of superheavy elements Z=119 
and Z= 120 + are calculated using ab initio approach. Core polarization and dominating correlation 
effects are included to all orders. Breit and quantum electrodynamic effects are also considered. 
Similar calculations for Cs, Fr, Ba + and Ra + are used to control the accuracy. The dependence of 
the hyperfine structure constants on nuclear radius is discussed. 
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I. INTRODUCTION 

Study of the hyperfine structure of heavy and super- 
heavy elements is an important source of the informa- 
tion about nuclear structure of these elements (see, e.g. 
Ref. [H, 0]). Hyperfine structure intervals are propor- 
tional to nuclear moments, such as magnetic dipole mo- 
ment, electric quadrupole moment, etc. The values of 
these moments can be extracted from the comparison of 
the calculations and the measurements. Apart from that, 
the hyperfine structure intervals are sensitive to electric 
charge and magnetic moment distributions within the nu- 
cleus. Parameters of these distributions can often be ex- 
tracted from the analysis of the hyperfine structure sub- 
ject to sufficient experimental data and the accuracy of 
the calculations. 

The hyperfine structure analysis can be even more im- 
portant for the superheavy elements [Z > 100) where 
sources of the information are very limited. The study of 
the superheavy elements are motivated by the hypothet- 
ical island of stability in the region Z=114 to Z = 126 
where shell closures are predicted (see, e.g., [3j). Ele- 
ments up to Z = 118, excluding Z = 117, have been 
synthesized (see, e.g., Refs. 0, 0j) and evidence for the 
naturally occurring element Z = 122 has been reported 

The use of the hyperfine structure analysis is limited 
to nuclei with odd number of protons or neutrons. The 
heaviest examples of such nuclei which can be found 
in the literature as being already observed include, e.g. 
fffUup and ffiUuq 0. There are numerous similar ex- 
amples for smaller Z. The information about nuclear 
magnetic dipole and electric quadrupole moments, charge 
and magnetic moment distribution for these elements is 
practically absent. 

Some standard approaches to the analysis of the hyper- 
fine structure do not work for very high Z . Consider, for 
example the Fermi-Segre Q formula with the Casimir 
relativistic factor. It expresses the hyperfine structure 
constant of the s-state of an external electron via its wave 
function in the origin. 

A s = const x \^(Q)\ 2 F re i(Za)(l - S)Ry, (1) 

where a is the fine structure constant, Ry is Rydberg, S 
is the correction due to finite nuclear size and F re \(Za) 



is the relativistic factor 

F re i(Za)= I 7 = V / 1 _ ~ZW. (2) 

7(47 - 1) 

Formulae I|ll2p are widely used in the hfs analysis, how- 
ever they fail at very high Z. It is easy to see that 
the relativistic factor ^ turns to infinity at Z = 118.7. 
Therefore, it is likely to overestimate the relativistic cor- 
rections at smaller Z as well. The reason for this is that 
the formulae treat finite nuclear size as a small correc- 
tion. Hydrogen-like wave functions for a point-like nu- 
cleus are used to calculate relativistic factor. However, 
it is known that the finite nuclear size correction for su- 
perheavy elements is not small and cannot be treated as 
a perturbation (see, e.g. 

A combination of an analytical and numerical ap- 
proaches were recently used in Ref. 9] to study the de- 
pendence of the hyperfine structure intervals on nuclear 
radius. A formula was suggested which is in good agree- 
ment with accurate numerical calculations for s-states of 
atoms wit Z < 100. However, this formula also fails at 
higher Z. 

In this paper we don't use any analytical approaches 
but just perform accurate numerical calculations of the 
hyperfine structure constants for superheavy atoms. We 
demonstrate that the calculation which use finite-size 
nuclei with realistic charge and magnetic moment dis- 
tribution are very similar to the calculations for lighter 
atoms. We consider elements El 19 and E120 + . The 
latter may have hyperfine structure if there is an iso- 
tope with odd number of neutrons. Neither of these ele- 
ments have been synthesized yet. However, the ways of 
their production and physics of their nuclei are discussed 
in literature [l(|. These elements are heavier than any 
known element. Therefore, if calculation of the hyper- 
fine structure brings no surprises for them one can ex- 
pect no surprises for lighter elements as well. Also these 
elements have very simple electron structure with one 
external electron above closed shells. Therefore, very ac- 
curate calculations are possible for the elements. In our 
previous work [Tl| we have calculated energy levels of 
E119 and E120 + . Apart from some expected relativistic 
effects like larger fine structure and stronger attraction 
of the s-states to atomic core, the spectra of these su- 
perheavy elements is very similar to the spectra of their 
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lighter analogies, Fr, Ra + , Cs and Ba + . We expect simi- 
lar trend for the hyperfine structure and we perform the 
calculations for the same set of atoms. This gives us an 
estimate of the accuracy of the results for superheavy el- 
ements. We stress that although the calculations bring 
no major surprise, the dependence of the hyperfine struc- 
ture on the nuclear radius is significantly stronger for the 
superheavy elements than for their lighter analogies. 



II. METHOD OF CALCULATION 

We perform the calculations using a totally ah initio 
method developed in our previous works 111, E EU E3j 
El, EH- K starts from the relativistic Hartree-Fork (RHF) 
calculations for atomic core and includes dominating cor- 
relation and all core polarization corrections to all orders. 

Single-electron orbitals are found by solving a system 
of the RHF equations for N — 1 electrons of the closed- 
shell core (the V N ~ 1 approximation). The RHF Hamil- 
tonian has a form 



H = cot • p + (p - l)mc 2 + KucO) + V 



N-l 



(3) 



Here a and (3 are Dirac matrixes, V nuc (r) is nuclear po- 
tential, y^ -1 = Vdir + Vcxch is the sum of the direct and 
exchange Hartree-Fock potentials, N is the number of 
electrons. At distances much larger than nuclear radius 
tat nuclear potential is given by V nuc {r) — —Ze 2 /r, at 
short distances V nuc (r) is obtained by numerical integra- 
tion of the Fermi distribution of nuclear electric charge. 
We use d = 2.3 fm as the^ thickness of the distribution 
and the data from Ref. [l7[ for the radii (see Table HJ. 
We use tat = 1.1(2.5Z) 1 / 3 fm for superheavy elements. 

The hyperfine interaction (HFI) is included in a self- 
consistent way as well. The time-dependent Hartree Fock 
method (TDHF) [12{ |. which is equivalent to the well- 
known random-phase approximation (RPA) is used for 
this. To take into account finite nuclear size we use a sim- 
ple model which represents the nucleus as an uniformly 
magnetized ball. In our calculations magnetic nuclear ra- 
dius is the same as the electric one. However, these two 
parameters can be varied independently. 

The HFI Hamiltonian is given by 



Hhfi = e/x-F(r), 



F(r) 



rxot 
— s — 



r < r n 
t^, r> r v 



(4) 



(5) 



where r m is the magnetic nuclear radius. 
The TDHF equations have a form 



(H - e a )8^ a = (-F z - 5V 1 "- 1 + 8e a )^ a 



Se a = {i) a \F z +5V N -^ a ) 



(0) 
(7) 



Here index a numerates states in the closed-shell core. 
These equations are solved self-consistently for all states 
in the core. 

States of the valence electron are calculated in the 
frozen field of atomic core complemented by the corre- 
lation potential operator S [l3[ 



(Ho + E - e)^° = 



(8) 



Here index v numerates valence states. Correlation po- 
tential £ includes all lowest second-order correlation cor- 
rections and dominating higher-order correlation correc- 
tions 0, EH- These higher-order correlations include 
screening of Coulomb interaction and hole-particle inter- 
action. They are taken into account in all orders. Solv- 
ing equations ijHJ) for valence states we find the so called 
Brueckner orbitals for the valence states. This is empha- 
sized by using superscript BO for the orbitals. 

Total energy shift for the valence state v due to HFI 
and correlations is given by 



Se v = ^ u \F z + SV^+St:\^ )- 



(9) 



Here 5T, is the change to the correlation potential S due 
to the hyperfine interaction. The term with 5T, is often 
called structure radiation. Finally, there is a contribution 
due to the renormalization of the many-electron wave 
function (see, e.g. [T^ ]) 

Se noTm = -(ip v \F z +5V N - 1 \ip v )(il> v \d£/dE\ip v ). (10) 

Magnetic dipole hyperfine structure constant A v for 
the valence state v is given by 



A v — 



fie 



Se v 



2m p I y/j v (jv + l)(2jv + l) 



(11) 



A. Breit and QED corrections 

It is hard to claim high accuracy of calculations for su- 
perheavy elements without considering Breit and quan- 
tum electrodynamics (QED) corrections. We include 
Breit corrections in a v ery accurate way described in 
our previous works [lfl [l8j]. The QED correction are 
included approximately via the QED potential suggested 
in Ref. [H|. 

The Breit operator has the form: 



h B = - 



oli ■ ol 2 + (oti ■ n)(a 2 • n) 
Yr 



(12) 



where r = nr, r is the distance between electrons, and 
a is the Dirac matrix. It correspond to the zero energy 
transfer approximation and includes magnetic interaction 
and retardation. 

Similar to the hyperfine interaction, Breit operator 
induces a correction to the self-consistent Hartree-Fock 
potential, which is taken into account in all orders in 
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TABLE I: Isotopes of Cs, Ba, Fr and Ra for which hyper- 
fine structure constants have been calculated in present work. 
Magnetic moments are in nuclear magnetons; gi = fj,/I. 



Isotope 




I 


9i 


r at [fin] 


rms[fm] a 


133 Cs 


2.582024 


7/2 


0.737721 


5.671 


4.8041 


135 Ba 


0.837943 


3/2 


0.558629 


5.703 


4.8273 


211 Fr 


4.00(8) 


9/2 


0.889 


6.717 


5.5545 


225 Ra 


-0.7348(15) 


1/2 


-1.4696 


6.887 


5.6781 



"Reference. |17| 



Coulomb interaction by iterating the RHF equations with 
the potential 



V 



N-l 



V c + V 1 



(13) 



where V c is the Coulomb potential, V B is the Breit po- 
tential. The same potential (fTB")) goes to the left and 
right-hand side of the TDHF equations © ■ 

For the QED corrections we use a radiative potential 
derived in Ref . [l9[ . This potential was chosen to fit accu- 
rate calculations of the QED corrections to the energies. 
It may give less accurate results for the hyperfine struc- 
ture. Therefore, we consider current calculations of the 
QED corrections as rough estimations only. 



III. RESULTS 

Table Q] lists isotopes of lighter analogies of the super- 
heavy elements E119 and E120 + for which the hyper- 
fine structure constants are calculated. The results of 
the calculations are presented in Table HU Here RHF 
corresponds to the (ip v \F z \ip v ) matrix elements with the 
Hartree-Fock wave functions ip v , the RPA corresponds 
to the (ip v \F z + 8V N ~ 1 \4) V ) matrix elements; BO and 
RPA(BO) columns correspond to the same matrix el- 
ements but with Hartree-Fock wave function replaced 
by Brueckner orbitals; the "Str+Norm" column includes 
structure radiation and renormalization. 

As can be seen from the table the most important cor- 
rections are the many-body corrections associated with 
the core polarization effect (RPA) and with the corre- 
lation interaction of the external electron with the core 
(BO). These effects follow approximately the same pat- 
tern when moving from light to heavy atoms. This means 
that the accuracy of the results should be about the same 
for all atoms and ions. 

Breit contribution is small and can be neglected in all 
cases. This is because Breit contributions are propor- 
tional to lower powers of Z than other relativistic effects. 
The QED corrections are large for s states. They reduce 
the hfs constants of these states by about 1%. 



We also study the dependence of the hyperfine struc- 
ture constants on nuclear radius. This is done numeri- 
cally by calculating the hfs constants at different radius 
and then calculating the derivative dA/dr n numerically. 
It is convenient to represent the results in a form of the 
dimensionless constant n^r as in Ref. Q 



Khr 



SA V / A v 

Sr n /r n 



(14) 



Here Ay is the hyperfine structure constant of the valence 
state v, r n is nuclear radius. The results are presented in 
Table Hill There are few things to note here. First, the 
effect in superheavy elements is much larger than in their 
lighter analogies. Second, the effect for s and pi/ 2 states 
is significantly different. This represents an opportunity 
to use the measurements of the hyperfine structure in su- 
perheavy elements not only to extract nuclear magnetic 
moments but also to get some information about nuclear 
radius. Note finally that the analytical formulae describ- 
ing the dependence of the hyperfine structure on nuclear 
radius presented in Refs. [3, Q do not work here. They 
are not just inaccurate, they give absolutely meaningless 
results. The reason for this is that the effect is large and 
cannot be treated perturbatevely. 

Calculating Khr we assume that magnetic and electric 
radii of the nucleus are the same. However, the pro- 
gram allows to treat them independently and calculate 
two partial derivatives dA/dr n and dA/dr m , where r n 
is electric radius and r m is magnetic radius. Such calcu- 
lations show that the hyperfine structure constants are 
more sensitive to the change of the electric radius. Cor- 
responding partial derivative is approximately two time 
larger that those over the magnetic radius. This is true 
for both s and P1/2 states. 



IV. CONCLUSION 

The hyperfine structure of lowest s and px/2 states of 
the superheavy elements Z=119 and Z=120 + have been 
calculated with an uncertainty of a few percent. The de- 
pendence of the constants on nuclear radius is presented. 
The results may be used for experimental studies of nu- 
clear, spectroscopic and chemical properties of the ele- 
ments. 
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TABLE II: Hyperfine structure constants of the lowest s 1 / 2 and P1/2 states of Cs, Fr, E119, Ba + , Ra + , and E120 + in different 
approximations in MHz (Cs,Fr,Ba+, Ra+) and gj x MHz (E119, Ef20+). 



Atom 


State 


RHF 


RPA 


BO 


RPA(BO) 


Breit 


Rad. 


Str+Norm 


Total 


Exp. 


Cs 


6s 


f425 


1718 


1970 


2325 


6 


-21 


-31 


2279 


2298.2 a 




6pi/2 


f6f 


202 


240 


294 








5 


299 


291.89° 


Fr 


7s 


5791 


6875 


7716 


8967 


33 


-162 


-120 


8718 


8713.9 6 




7P1/2 


623 


772 


968 


1180 


-4 


-4 


8 


1180 


1142 b 


E119 


8s 


39344 


46781 


45531 


53306 


210 


-553 


-1315 


51648 






°Fl/2 


5141 


6165 


8751 


10506 


-65 


-41 


135 


10535 




Ba+ 


6s 


2607 


3095 


3147 


3684 


8 


-42 


-82 


3568 


3591.6706(3) c 




6^1/2 


441 


530 


568 


674 


-1 





4 


677 


664.2(2.4)" 


Ra+ 


7s 


-21357 


-25022 


-25114 


-28986 


92 


436 


668 


-27790 


-27684(13) d 




7P1/2 


-3626 


-4330 


-4746 


-5611 


19 


19 


-31 


-5604 


-5446(7)"* 


E120+ 


8s 


74195 


86640 


80396 


92884 


352 


-837 


-2790 


89609 






8P1/2 


16883 


19849 


22286 


26218 


-117 


-92 


12 


26021 





"Reference. 2£| 

^Reference. 2U 

"Reference. 2^ 

d Reference. 23( 



TABLE III: Sensitivity of the hyperfine structure constants 
to the change of nuclear radius («fe r ). 



Cs 


Fr 


E119 


E20 4 




6s 


7s 


8s 8pi /2 


8s 


8pi/2 


-0.024 a 


-0.11" 


-0.45 -0.28 


-0.46 


-0.22 



"Reference. [9( 
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